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We review recent experiments on entanglement, Bell's in- 
equality, and decoherence-free subspaces in a quantum regis- 
ter of trapped 9 Be + ions. We have demonstrated entangle- 
ment of up to four ions |l| using the technique of M0lmer and 
S0rensen fj. This method produces the state (|TT) + I4I))/V2 
for two ions and the state (| UU) + | TTTT))/V2 for four 
ions. We generate the entanglement deterministically in each 
shot of the experiment. Measurements on the two-ion entan- 
gled state violates Bell's inequality at the 8a level ||. Be- 
cause of the high detector efficiency of our apparatus, this 
experiment closes the detector loophole for Bell's inequality 
measurements for the first time. This measurement is also 
the first violation of Bell's inequality by massive particles 
that does not implicitly assume results from quantum me- 
chanics. Finally, we have demonstrated reversible encoding 
of an arbitrary qubit, originally contained in one ion, into a 
decoherence-free subspace (DFS) of two ions Q|. The DFS- 
encoded qubit resists applied collective dephasing noise and 
retains coherence under ambient conditions 3.6 times longer 
than does an unencoded qubit. The encoding method, which 
uses single-ion gates and the two-ion entangling gate, demon- 
strates all the elements required for two-qubit universal quan- 
tum logic. 



I. INTRODUCTION 

Trapped ions are a promising candidate system for im- 
plementing quantum computing (QC), with potential ap- 
plications like efficient factorization of large numbers |]] 
and efficient searching of large databases Q . While large- 
scale QC is still far off, trapped-ion quantum registers 
have already demonstrated all the essential ingredients 
of QC, including initialization to a known quantum state 
|7|,p|p0|] , efficient detection of final states p|,^0| , long qubit 
coherence times jSO^Mj, an d a universal set of quan- 
tum logic gates fl] jl0| , |14| , [15f ■ Recent work on trapped-ion 
QC at NIST has used this basic set of tools to produce 
four-particle entanglement Jlj , observe violation of Bell's 
inequality || , and encode quantum information into a 
decoherence-free subspace j|. 

In the experiments described below, the qubits are 




FIG. 1. Partial level diagram of the 9 Be+ ion, showing the 
cycling transition for er_ resonant light and the two-photon 
off-resonant coupling for Raman transitions. Wavelengths for 
all transitions are approximately 313 nm. 



electronic states of 9 Be ions held in a linear RF trap 
of the type described in reference [[T^| . The relevant elec- 
tronic energy levels of the 9 Bc + ion are shown in Fig. 1. 
Each ion encodes one qubit in its 2 Si/ 2 \F = 2, nip = —2) 
and 2 Si/ 2 \F = l,mj? = —1) hyperfine sublevels, denoted 
by ||) and ||), respectively. The ions are readily initial- 
ized in ||) by optical pumping. We detect the number of 
ions in state ||) by applying left circularly polarized laser 
light to excite the ||) - 2 Ps/2 transition at about 313 nm. 
The polarization of the light ensures excitation to the 
2 P 3 / 2 |-F = 3,mp = —3) hyperfine sublevel, which can 
only decay to the ||) sublevel by dipole selection rules. 
Thus an ion in state ||) can scatter thousands of photons 
without changing its hyperfine state. On the other hand, 
state |t) is spectrally resolved from ||) by the hyperfine 
splitting of 1.25 GHz and therefore scatters very little 
light. The resulting detection accuracy per ion is about 
98% in a single shot of the experiment. 

We produce controlled rotations of a single qubit by 
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driving stimulated Raman transitions between states ||) 
and | |). The laser beams driving the transitions are 
tuned about midway between the 2 P3/2 and 2 Pi/2 states. 
If the frequency difference is^is between the Raman laser 
beams is equal to the energy splitting vq between states 
|t) and ||), we observe Rabi flopping between |f) and ||) 
corresponding to the evolution 

l4)-cos^||)+ e ^sin^|T) (I) 

IT) - cos t |t) -e-« sin 

Here is proportional to the Raman pulse duration and 
cj> is the phase difference between the Raman beams at 
the position of the ion, referred to as the "ion phase" . 
For two ions, we write the ion phases as fa, fa. We can 
effect common-mode changes of fa and fa by changing 
the phase of the RF synthesizer that controls the Ra- 
man difference frequency. By changing the strength of 
the trap, we can control the ion spacing precisely, allow- 
ing us to make differential changes to fa and fa. These 
two techniques combined allow us to control fa and fa 
independently for two ions, so that we can perform single- 
qubit rotations on one ion without affecting the other M . 
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FIG. 2. Level diagram for the S0rensen-M0lmer entangling 
gate in the case of two ions. Here, as opposed to Fig. 1, arrows 
indicate two-photon Raman transitions. Raman transitions 
near both red and blue sidebands are applied simultaneously 
to create the state (| j,J.) + | tt))/V2- 



II. GENERATING ENTANGLEMENT 



The potential created by the ion trap is harmonic with 
axial frequency vqm of about 5 MHz. Doppler cooling al- 
lows us to cool the ions to temperatures T ~ /wcm/^b- 
The Doppler-cooled ions are strongly coupled by the 
Coulomb interaction and form a rigid crystal, which in 
these experiments is a string of ions lying along the trap 
axis. Because the ions are strongly coupled, their quan- 
tized motion is best described in terms of normal modes 
with frequencies z^. For two ions, for instance, the mo- 
tion of the ions along the axis decomposes into symmetric 
and antisymmetric normal modes at frequencies ucm and 
VSvqm- The lowest- lying normal mode for any number 
of ions is the center-of-mass mode at frequency vqm- In 
this mode the ion string moves as a unit. 

We can control the motional quantum state of the ions 
by driving stimulated Raman transitions that couple spin 
and motion. In our geometry, the wavevector difference 
between Raman beams lies along the axis of the ion trap, 
so the Raman beams exert a dipole force only along the 
axis. Setting Vdis equal to vq + mv% (m an integer) in- 
duces Rabi flopping between states \l)\n) and |t)|n + m) 
with evolution like that of Eq. (Q). Here \n) denotes the 
nth Fock state of the ith normal mode. In particular, 
we can cool the ions by driving a pulse with 9 = n at 
i^diff = vq — Vi (the "red sideband" ) and optically pump- 
ing from 1 1) back to | |). Repeating the cycle several 
times permits preparation of the ion string in its mo- 
tional ground state 



We have entangled strings of up to four ions |f| us- 
ing the technique proposed by S0rensen and Maimer || . 
In this method, one drives two Raman transitions si- 
multaneously, with the same Rabi frequency on all ions. 
Fig. 2 shows an energy-level diagram including the rel- 
evant transitions. The detuning S ensures that popula- 
tion transfer between motional states is nonresonant and 
therefore small. However, the sum of the two Raman 
transition frequencies is 2u^ } so that for any pair of ions 
the process |||) «-» |tt) is resonant. For N ions initially 
in |t), the resulting evolution produces the state 

\fa) = (ltt) + ltt))/V2 _ (two ions) 

(|UU>+ |TTTT»/V2 (four ions) {! 

both of which are maximally entangled states. 

To implement this technique, we pass one of the Ra- 
man beams through an electro-optic modulator (EOM) 
to generate the two required Raman difference frequen- 
cies and uniformly illuminate the ion string with both 
Raman beams. For two ions, we tune the Raman dif- 
ference frequencies close to the motional sidebands of 
the antisymmetric mod e; for four ions, we use the mode 
with frequency y / 29/5 vcm- In both cases, the mode 
we choose has equal amplitudes of motion for all ions, 
so that the Raman Rabi frequencies for all ions are the 
same under uniform illumination. 
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In general, the state produced by this method is char- 
acterized by some density matrix p that approximates 
the desired state \ip e ). To show entanglement of the ex- 
perimental state p, it is sufficient to show that the fidelity 
F = (ipe\p\i/) e ) exceeds 1/2 fli. We have 



F = (pii-a + Pn-rr)/ 2 + Ipu-ttI 



(3) 



We measure the diagonal density matrix elements 
Pi— 1,1— J.! PT---T,T~ T by detecting the number of ions in 
\l). To find the far-off-diagonal coherence |px...x|...t|, we 
rotate all ions together as in Eq. (|l|) with 8 = ir/2 and 
with phases fa — <j> relative to the entangling pulse. Then 
we measure the parity II of the ion string, defined as +1 
for an even number of ions in ||), —1 for an odd number 
in ||). It is easy to show that 



n = 2|p u ... TT |cosiV</> 



(4) 



for N ions, if all other off-diagonal coherences are absent. 
Sweeping the phase of the RF synthesizer that generates 
the Raman difference frequency gives a sinusoidal signal 
of II vs. fa as shown in Fig. 3. We extract the coherence 
\Pll— Ttl ^ fitting a sine wave to the data of Fig. 3. The 
density matrix elements needed to calculate the fidelity 
are given in Table 1 for two and four ions. For two ions we 
find F = 0.83 ± 0.01, while for four ions F = 0.57 ± 0.02, 
demonstrating entanglement in both cases |Q. It is im- 
portant to note that the data used to calculate F are not 
post-selected. Since the evolution is deterministic, our 
method produces an entangled state in each shot of the 
experiment. 



TABLE I. Partial density matrix of the state produced by 
the entangling pulse. From Eq. (|^), only these matrix ele- 
ments are needed to determine whether the experimentally 
generated state is entangled. Errors are ±0.01 for two ions 
and ±0.02 for four ions. 
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FIG. 3. Measurements of the parity II as a function of RF 
synthesizer phase <f> for two and four ions. The data show clear 
sinusoidal dependence as in Eq. (^|) , allowing us to extract the 
far-off-diagonal coherence |pu---Ttl- 



assume local realism. Bell ]19| ] proved that quantum me- 
chanics is actually inconsistent with local realism. Quan- 
titative conditions for falsifying local realism by experi- 
mental measurements were first proposed by Bell jy| and 
were developed further by Clauser, Horne, Shimony, and 
Holt (CHSH) (^0|. In such an experiment, one prepares a 
pair of particles in a repeatable way, subjects them to in- 
dependent classical manipulations parametrized by vari- 
ables $i, $2, and measures a classical property yielding 
a value Pi($i), P2(3b) = ±1 for each particle. Defining 
the correlation function g(<E»i , <J? 2 ) = -Pi($i)-P2($2), one 
finds that the Bell inequality 

£(ai,<5i,/?2,72) = k(<5i,72) - 2(ai,72)| (5) 
+ \q(S 1 ,(3 2 )+q(a 1 ,p 2 )\<2 

holds for any values a\, 8% (resp. (3 2 , 12) of $1 (resp. $2) 
under the assumption of local realism. 



III. VIOLATION OF A BELL INEQUALITY 

Classical physics obeys the principle of local realism, 
which states that objects have definite properties whether 
or not they are measured and that measurements of an 
object's properties are not affected by events taking place 
sufficiently far away. As shown by Einstein, Podolsky, 
and Rosen p8[, quantum mechanics is incomplete if we 



In our experiment we apply the entangling oper- 
ation of Eq. (Q) to two ions to prepare our initial state. 
Then we perform a carrier rotation as in Eq. ([!]) with 
9 = tt/2 and phases fa, fa relative to the entangling pulse 
phases. Finally, we detect the number of bright ions. The 
calibration of fa and fa against the trap strength and the 
RF synthesizer phase were performed using assumptions 
from classical physics only ||. We can therefore iden- 
tify fa, fa with the classical control parameters <!>i,<I>2, 
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so the experiment is of the type proposed by Bell and 
CHSH. Quantum mechanics predicts that, for an ideal 
realization of our experiment, B(a.i, Si, 02, 72) attains a 
maximum value of 2\/2 for 



ai = -tt/8 <Si = 3tt/8 
P2 = -tt/8 72 = 3tt/8 



(0) 



contradicting the bound B(ot\, <5i, /3 2 , 72) < 2 required by 
local realism. 



We repeated the experiment TVtot = 20 000 times for 
each of the four sets of phases given by Eq. (|). These 
data are presented as histograms in Fig. 4. We see 
that the three cases of zero, one, and two ions bright are 
readily distinguishable. The overlap of histograms corre- 
sponding to the three cases gives an inaccuracy of 2% in 
discriminating between the three cases. This inaccuracy 
mostly arises from pumping of the | J.) state to the 1 1) 
state because of imperfect polarization of the detection 
light. Counting the number of events Nq (Ni,N2) with 
zero (one, two) ions bright, we extract the correlation 
function using 



{Nq + N 2 ) - Ni 

JVtot 



(7) 



and obtain a set of four values of q((f>i, <fi 2 )- We repeated 
the entire procedure five times to ensure that the mea- 
sured Bell signal did not vary significantly from run to 
run. The five sets of correlation functions and the re- 
sulting Bell signals are shown in Table 2. The statistical 
error of each Bell signal was ±0.01. Averaging the Bell 
signals together, we find 



B(- 



7T 3-7T 



3tt 



) = 2.25 ±0.03 



(8) 



an 8(7 violation of the Bell inequality Eq. (|5|). The error 
in the average Bell signal significantly exceeds the error 
of each individual measurement of the Bell signal. This 
discrepancy arises from errors in setting the phases of 
Eq.®. 




FIG. 4. Histograms of fluorescence signals for the four sets 
of phases given in Eq. (^j). The vertical arrows break up the 
data into cases of zero, one, or two ions bright with 98% 
accuracy. 



TABLE II. The correlation values and resulting Bell's sig- 
nals for five runs of the experiment. The phase angles 
Qi, Si, P2, 72 are given by Eq. (|^). The statistical errors are 
0.006 and 0.012 for the values of q and B respectively. 



q{ai,(5 2 ) g(ai,7 2 ) q{8i,f3 2 ) 9(^1,72) B(a ly 61, (3 2 , 72) 
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2.275 
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0.559 


-0.551 
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0.541 


0.596 


0.537 


-0.571 


2.245 



each other's lightcones, in principle the detections could 
influence each other, leading to spurious correlations. 



IV. ENCODING INTO A DECOHERENCE-FREE 
SUBSPACE 



This experiment gives the first violation of the Bell in- 
equality for massive particles obtained using a complete 
set of correlation measurements. A previous experiment 
pl[ deduced a violation from an incomplete set of mea- 
surements. Furthermore, our data uses the outcome of 
every shot of the experiment, so the violation is obtained 
without the use of the fair-sampling hypothesis. In prin- 
ciple, sampling or post-selection of the data allows vio- 
lation of Eq. (§) even by local hidden-variable theories, 
the so-called "detector loophole" P4] . Our experiment 
closes the detector loophole for the first time. However, 
the "locality loophole" |23| remains open for our data; 
since the detection events on the two ions occur within 



Information stored in a quantum system is liable to be 
lost through coupling of the system to its environment, 
a process called decoherence |24|,^5]]. In particular, many 
proposed quantum computer implementations decohere 
through a mechanism that couples each qubit to the envi- 
ronment equally strongly p6| , p7f . Encoding quantum in- 
formation into a decoherence-free subspace (DFS) p^-|3C|] 
whose states are invariant under coupling to the environ- 
ment allows information storage even in the presence of 
this type of decoherence. We have demonstrated the im- 
munity of a DFS of two atoms to environment-induced 
dephasing and implemented a technique for encoding an 
arbitrary physical qubit state into the DFS. 
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The type of decoherence considered here is collective 
dephasing, in which each qubit undergoes the transfor- 
mation | J.) — > ||), |t) — * e 4 ''!?) f° r C an unknown phase 
that is the same from qubit to qubit. One can construct 
a DFS of two qubits that protects against collective de- 
phasing from the basis states |||), |IT)- Any superposi- 
tion of these states is clearly invariant under collective 
dephasing, and since there are two states, we can encode 
one "logical qubit" in this subspace of the two "physical 
qubits" . To demonstrate our encoding technique, we first 
prepare a state of form 

\l)(a\l) + 6|T)) M 2 + |&| 2 = 1 a, b complex (9) 

where one of the two physical qubits initially contains 
the quantum information. In general, we can prepare 
this state by driving the carrier transition of Eq. (|l|) on 
both ions simultaneously, once with 9 = /3, and again 
with 9 = /3 and <fri shifted by tt. We set 4>2 = ol for both 
pulses. The final state has a = cos 2/3, b = e la sin 2/3. 
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FIG. 5. Decay of DFS-encoded state (circles) and test state 
(crosses) under engineered dephasing noise. The noise is ap- 
plied for a fraction of the delay time of about 25 fis between 
encoding and decoding. Coherence data are normalized to 
their values for zero applied noise. The fit lines are exponen- 
tial decay curves. The DFS data is fit for comparison. 



We encode the state of Eq. (g) into the DFS in two 
steps. First we apply the inverse of the entangling gate, 
yielding a(|||) - «|TT» + HUT) - i\U))- Then we drive 
the carrier with 9 = ir/4,<pi = n/2,(j)2 = to obtain 
IV'dfs) = o(UT) + *|U» + HUT) - i|U»- The informa- 
tion stored in the physical qubit of ion 2 is now encoded 
in the logical qubit of the DFS. In the experiment we 
take \a\ — \b\, though our method permits preparation 
and encoding of any state of the more general form. 



To read out the encoded information, we reverse the 
carrier pulse in the encoding and apply the entangling 
gate to decode IV'dfs) into | j) + b|t))- After de- 
coding, we rotate ion 2 as in preparing the state of Eq. 
(^) but with the phase on ion 2 changed to a' . We then 
measure the probability P2 of finding both ions in ||). 
P2 varies sinusoidally with a — a' , and the magnitude of 
oscillation is equal to the coherence, C, of ion 2. Since 
we set |a| = |6| = 1/V2, ideally C = 1. Departures from 
C = 1 measure the effects of both decoherence and im- 
perfect logic. We verified that C is independent of a, 
thus showing that our encoding method works even if we 
have no information about the phase of the input state. 



To study the effects of decoherence on the DFS- 
encoded state, we engineered a collective dephasing en- 
vironment using an off-resonant laser beam with a ran- 
domly varying intensity. The beam induces shifts of vq 
that are common to both ions through the AC Stark ef- 
fect. The DFS state should resist the dephasing effect 
of this environment. The coherence of ion 2 in the test 
state ||) (||) + e l< ^\\))/y/2, however, should be sensitive to 



collective dephasing. We measure the decay of the test 
state by simply turning off the encoding and decoding 
sequences in the procedure used to measure the decay of 
the DFS-encoded state. We applied the noise beam to 
the test state and encoded state during a delay time of 
about 25 /us, as shown in Fig. 5. The coherence without 
applied noise is about 0.69 for the test state and about 
0.43 for the encoded state; they depart from 1 because of 
imperfect logic gates and detection. For white-noise in- 
tensity fluctuations of the Stark-shifting beam, we expect 
C to decay exponentially for the test state, as shown by 
the fit line. We also fit the coherence of the DFS state to 
an exponential decay for comparison. The decay rate of 
the test state is 0.18(1) /is -1 , while the decay rate of the 
DFS state is 0.0035(50) /is -1 , consistent with zero decay. 

We have also measured the storage times of the en- 
coded and test states under ambient conditions in our 
laboratory. The data are shown in Fig. 6. Here we 
leave a variable delay time between encoding and decod- 
ing to give the ambient noise time to act. The decoher- 
ence of the test state is dominated by ambient fluctu- 
ating magnetic fields whose frequencies lie primarily at 
60 Hz and its harmonics. These fields randomly shift loq 
through the Zeeman effect. Since these fields are roughly 
uniform across the ion string, they induce collective de- 
phasing similar to that created by the engineered en- 
vironment. We empirically find the decay of both test 
and encoded states to be roughly exponential, as indi- 
cated by the fit lines. The decay rate of the test state 
is 7.9(1.5) x 10~ 3 /is^ 1 , while the decay rate of the DFS 
state is 2.2(0.3) x 10 _3 jUS -1 . The long lifetime of the DFS 
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FIG. 6. Decay of DFS-encoded state (circles) and test state 
(crosses) under ambient decoherence. We vary the delay time 
between encoding and decoding to give the ambient noise a 
variable time to act. Coherence data are normalized to their 
values for zero delay time. The fit lines are exponential de- 
cay curves for purposes of comparison and are not theoretical 
predictions. 



state shows that collective dephasing from magnetic field 
noise is the major ambient source of decoherence for the 
test state. The loss of coherence of the encoded state is 
consistent with degradation of the decoding pulses [p|j3l| 
due to heating of the ion motional state over the delay 
time Il7 . 
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FIG. 7. Schematic of the proposed "quantum CCD" ar- 
chitecture for a large-scale ion-trap quantum computer. The 
"logic ions" , which encode the quantum information, are held 
in a storage region. To perform a logic operation, we shuttle 
the relevant logic ions into the "accumulator" region, where 
they interact with lasers. Auxiliary "cooling ions" are pro- 
vided for sympathetic cooling of the logic ions, which will 
almost certainly heat up during transfer to the accumulator. 



V. CONCLUSION AND OUTLOOK 

We have demonstrated four-particle entanglement, 
violation of Bell's inequality, and encoding into a 
decoherence-free subspace using a small quantum reg- 
ister of trapped ions. While these experiments are im- 
portant steps on the road to quantum computation, the 
major challenge remains: scaling up the techniques de- 
veloped for a few ions to build a large-scale quantum 
computer with thousands or millions of qubits. Fortu- 
nately, we have a plausible roadmap to our goal, the 
"quantum CCD" model of ion-trap quantum computing 
(Fig. 7). We envision a large ion trap with many in- 
dependent DC electrodes, which could be used to move 
ions between storage regions and "accumulator" regions 
in which quantum logic operations take place. Such a 
trap could hold thousands of ions, but avoids the slow 
gate speeds and detection inefficiency associated with 
large ion crystals by allowing us to work with only a few 
qubits at a time. Though the ions would most likely heat 
up while being moved around, the addition of another ion 
species would enable us to sympathetically cool the ions 



before performing logic |32f| . The DFS encoding demon- 
strated above would permit us to move qubits around 
without incurring dephasing from the spatial variation 
of magnetic fields across the device. In pursuit of the 
"quantum CCD" goal, we have built a trap supporting a 
double-well potential with spacing about 1 mm and plan 
to explore ways of maintaining quantum coherence and 
entanglement over these large distances. 
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